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1. Introduction 

The chiral quark condensate (QC) (yy) plays a central role in the non-perturbative sector of 
QCD. It is expected to be the order parameter which controls the spontaneous chiral symmetry 
breaking and sets the scale of the pion masses, provided its value is non-vanishing and its order 
of magnitude appropriate. To verify these basic requirements, high precision estimates of the QC, 
though welcome, are not necessary. 

At the same time there is no easy experimental access to the QC, as some theoretical assump- 
tions always enter. This suggests to look for different methods based on independent hypothesis 
which, to some extent and indirectly, can be tested in this way. Two examples are the classical 
measures of the QC extracted from the pion scattering length through chiral perturbation theory 
[jl|] or from the nucleon and B* — B mass-splitting sum rules 

On the other hand lattice QCD can give an estimate of the QC based on first principles. The 
simplest approach to compute the QC on the lattice uses the GMOR formula, which is based on the 
axial chiral Ward identity [Q]. Alternative determinations of the quark condensate on the lattice 
have been obtained in the framework of the £-expansion of QCD in a small volume and from 
the study of the Goldstone pole contribution to the pseudoscalar quark Green function 

The method proposed in [[J and summarized in this talk starts from the operator product 
expansion of the quark propagator at short euclidean distance: 

S(* 2 ) = ^(x 2 ) + ^ 2 C m {x 2 ) 5-JL C w (x 2 ) (W) + ■■■ (1.D 

The quark propagator on the l.h.s. is computed on the lattice and the QC is extracted fitting the 
r.h.s. in the chiral limit. Our final result reads (v>n//) MS (2GeV) = -(265±23MeV) 3 . 

The idea is simple and provides a non-perturbative test of the hypothesis behind, namely the 
validity of the operator product expansion in QCD. The precision achieved in this exploratory 
study is not high, but can be improved in time, being mainly computer-limited. However the 
method itself is flexible and powerful: in priciple, from the fit of eq. (pTT|), the quark masses and 
the gluon condensate which appears at G(x 2 ) can be obtained too. Another advantage is that the 
renormalization procedure is rather simple: once the quark propagator is renormalized, its operator 
product expansion is expressed in terms of renormalized quantities too. The same holds for the 
0(a) improvement. 



2. Lattice simulation 

We generated 180 gauge configurations in the quenched approximation with the f?(a)-improved 
Wison action, on a volume 32 3 x70 at j8 = 6.45. This value of the coupling corresponds to an in- 
verse lattice spacing aT x = 3.87(19)GeV, evaluated from the K and K* meson masses with the 
lattice plane method. The simulated values of the hopping parameter are K = 0.1349, 0.1351, 
0.1352, 0.1353 and correspond to light quark masses in the range m s /2 ^ m < m s . The values of 
the renormalized masses used in the chiral extrapolations and in eq. (1.1) have been taken from []7|]. 



The statistical errors are computed with the jackknife technique. The quantities directly evaluated 
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Figure 1: The bare form factors T.\ (left panels) and T.2 (right panels). From top to bottom: form factors 
in the interacting theory at k=0.1349; in the free lattice theory at infinite volume and in the chiral limit; 
corrected form factors (as defined in the text). 



from the numerical simulation are the scalar form factors of the quark propagator, defined from 

S(x 2 ) = j^Mx 2 ) + ^ 2 (x 2 )- (2-1) 

We have calculated the Wilson coefficients Cj(x 2 ) in eq. ( JTTT] ) at &{a s ) and resummed at the NLO 
in the MS scheme. The NLO expressions of the coefficients and of the evolution functions can be 
found in [g]. From here on, x indicates the euclidean distance in lattice units squared, (x/a) 2 . 

The scalar form factors Zi^^ju) are calculated in the Landau gauge and shown in fig. |] 
(top) for K = 0.1349 as function of x. The points show fish-bone curves due to the spread of the 
results obtained from different lattice sites which correspond to the same x 2 in the continuum limit, 
especially at short distance. Similar patterns can be observed in the form factors computed in the 
free lattice theory and plotted in fig. |I](center). This suggests we are observing lattice artifacts. We 
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can reduce the discretization errors from 0{a ) to &{a s a 2 ) by correcting the lattice data with the 
lattice artifacts evaluated in the free theory. In detail, we define 

(lf ont (x) \ 1 

2r r M = ^M-4"™,M ^"«W = o. (2-3) 

The result of the correction is shown in fig. [j](bottom). These are the data actually used in the 
extraction of the QC. 

3. Renormalization 

As anticipated, we need to compute only the renormalization constant of the quark field. For 
this purpose, we use the non-perturbative renormalization scheme named X-scheme in [||]- [||] and 
defined by the condition 

Z*(ji = l/x)Z 1 (x)=Vff ne (x) (3.1) 
in the Landau gauge and in the chiral limit. From eqs.(|l.l[) and ( |2.1| ), 

El(*) = 2^C/(*) + -. ( 3 - 2 ) 

In order to reach the chiral limit, linear and quadratic chiral extrapolations have been considered 
and the differences are included in the systematics. The result for Z^(/x = l/x) as obtained from 



eq. (3.1) is shown in fig. The final estimate of the renormalization constant is obtained from a 



constant fit in the range [9,25]: 

Z*(2GeV) = 0.871 ± 0.003 rtflf ± 0.020^ (3.3) 

This is also the of value of Z^ s (2GeV), in the Landau gauge. A fit which includes a long distance 
region, x E [9,40], has been also performed in order to evaluate the systematics. 

4. Determination of the quark condensate 



From eqs.(|l.l|) and (gjj) 



^(^M) = ^Cr S (x,M)m MS ( Ai )- i ^C^(x, Ai )(w) MS (M)^ 2 + --- (4-1) 

Therefore in the chiral limit the QC represents the leading term in the OPE of ^C^aO- A linear 
and a quadratic extrapolations have been performed in order to evaluate the systematics. 

In order to identify directly the QC as the intercept, the fits in x are performed on the normal- 
ized quantity 

~^0m -™ w + ^ 2> (4 - 2) 

Some results are collected in table [j] and the result of a constant fit in [9,25] is shown in fig. ||. 
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Figure 2: Values ofZy(fi = 2GeV) for different values of x = (x/a) 2 . The solid lines indicate the results 
obtained from a constant fit in x and the estimated systematic error. The dashed vertical lines show the x- 
range where the constant fit is performed. The result of a linear fit is also shown (dashed line); the estimate 
ofZy(fl = 2GeV) is given by the intercept. 



Constant fit [9,25] 


Linear fit [9,40] 


-(265 ±5) 3 
-(266 ±4) 3 


-(265 ±7) 3 



Table 1: Values of the chiral quark condensate in the MS scheme at the scale /i = 2GeV. 



As a check of the result we also extract the QC from E2(x,/l) by inverting the order of the 
fits: first, the fit in x including the mass term, then the chiral fit. We obtain a perfect agreement 
of the results with the central values given in table |T[ We note, however, that in order to stabilize 
the results of the fit, we have not been able to treat the quark masses as free parameters, but their 
values have been kept fixed to the ones determined in ref.[|7|]. Averaging the results and adding the 
systematic error we obtain 

(Wp(2GeV) = - (265 ± 5^ ± 22. vy ,. f MeV) 3 . (4.3) 

The systematic uncertainty amounts to 26%, the main contribution being the spread of the points 
in the fitting regions (18%). The reason is that the fitting intervals have been chosen in the region 
x G [9,40] in order to satisfy the condition 

1<v^^1/(«Aqcd) *=Q 2 (44) 

where the discretization errors are under control (lower limit) and the perturbative calculation of the 
coefficients at a scale pL ~ l/x is reliable (upper limit). But at the value of the lattice spacing con- 
sidered in the present simulation, this region is small and affected by non-negligible discretization 
effects, even after the correction of the lattice artifacts described in sec|2| has been applied. 

We estimate that neglecting higher order terms in the OPE introduce an uncertainty of ~5%. 
All the other sources of systematic error are due to the lattice uncertainties: 15% the error due 
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Figure 3: Values of 



'x,li)x 2 /4N C 



(\jf\jf) (p) + ff(x) as a function of x = (x/a) . A constant and a 



linear fit in x are shown. The absolute value of the quark condensate is given by the intercept. 



to the determination of the lattice spacing, 10% the quenched estimate of Aqcd, 9% the use of 
linear/quadratic chiral extrapolations, 5% the choice of the fitting interval, < 1% the error on renor- 
malization constant of the quark field and the use of vector or axial definitions of the quark masses. 
We could not quantify the finite volume effects, which should be small in the fitting regions (see 
[||] for a study in the free theory), and, above all, the quenching approximation. 
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